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ABSTRACT 

We present a unified model for the structure and appearance of accretion 
powered sources across their entire luminosity range from galactic X-ray binaries 
to luminous quasars, with emphasis on AGN and their phenomenology. Central 
to this model is the notion of MHD winds launched from the accretion disks 
that power these objects. These winds provide the matter that manifests as 
blueshifted absorption features in the UV and X-ray spectra of a large fraction 
of these sources; furthermore, their density distribution in the poloidal plane 
determines the "appearance" (i.e. the column and velocity structure of these 
absorption features) as a function of the observer inclination angle. This work 
focuses on just the broadest characteristics of these objects; nonetheless, it pro- 
vides scaling laws that allow one to reproduce within this model the properties of 
objects spanning a very wide luminosity range and viewed at different inclination 
angles, and trace them to a common underlying dynamical structure. Its general 
conclusion is that the AGN phenomenology can be accounted for in terms of 
three parameters: The wind mass flux in units of the Eddington value, m, the 
observer's inclination angle 6 and the logarithmic slope between the O/UV and 
X-ray fluxes ctox- However, because of a significant correlation between otox and 
UV luminosity, we conclude that the AGN structure depends on only two param- 
eters. Interestingly, the correlations implied by this model appear to extend to 
and consistent with the characteristics of galactic X-ray sources, suggesting the 
presence of a truly unified underlying structure for accretion powered sources. 
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1. Introduction 

The notion of AGN as an astronomical object of solar system dimensions and luminosity 
surpassing that of a galaxy has been with us for about half a century now. Since then, the 
advent of novel observational techniques, the accumulation of data and theoretical modeling 
has refined and advanced our notions as to what constitutes an AGN, with accretion onto 
a black hole as the source of the observed radiation now being universally accepted. At 
the same time, the discovery of galactic bright X-ray binary (XRB) sources, powered also 
by accretion onto compact objects (neutron stars and stellar size black holes) has extended 
the notion of accretion powered source to the stellar domain. Indeed, the general similarity 
of the X-ray spectral properties of AGN and galactic black hole can didates (GB HC) and 



XRBs in general, including their broad Fe Ka fluorescence features (IMillerl 120071 ) . argues 



for near horizon structures which are very similar, despite the hu ge disparity in the objects 



scales. This structure is thought to consist of a Shakura-Sunyaev ( iShakura fe Sunyaevlll973l ) 
disk that extends to the ISCO (innermost stable circular orbit) of the corresponding flow, 
supplemented by an overlying hot, X-ray emitting corona. 

Even though it is generally accepted that the AGN radiant energy is released by the 
accretion of matter onto a black hole (or in certain cases by extraction of the hole's rotational 
energy) in a region comparable to its horizon, there is plenty of evidence that a significant 
fraction of the AGN power is emitted, after reprocessing, at much larger radii. [One should 
note however, that accretion energy can also be transported outward not only radiatively 
but also me chanically by the viscous str esses that transport the accretion flow's angular 
momentum ( Blandford fe Begelman 19991 )]. Thus, the UV and optical lines that constitute, 



typically, a fraction / ~ 10% of the AGN bolometric luminosity, are emitted presumably 
by clouds at distances ~ 0.1 — 10 pc that cover a fraction / of the AGN solid angle. In 
addition to the line emission, the AGN ionizing continuum is also reprocessed into IR and 
far-IR radiation by matter at even larger distances, which apparently subtends an even larger 
fraction of the AGN solid angle (~ 50%). The geometry of this component is thought to be 
cylindrical (rather than spherical) with a column density that depends strongly on the angle 
9 of the observers' line of sight (LoS) with the symmetry axis. It was proposed that such 



a geo metry nicely unifies the Seyfert-1 and Seyfert-2 AGN subclasses (lAntonucci fe Miller 



1985) and also those of the broad and narrow line radio galaxies (BLRG - NLRG) (IBarthel 



19891 ). according to the angle 9: Thus, Seyfert-ls (or BLRG) are AGN in which the observer's 



LoS makes a small angle with their axis of symmetry, the column of the intervening cold 
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gas is small (Nh < 10 21 cm -2 ) and the continuum source and its surrounding broad line 
emission (concentrated in the inner AGN regions) are directly visible. Seyfert-2s (or NLRG) 
on the other hand, represent the same objects viewed at a large inclination angle, along 
which the column density to the source is much larger (Nh > 10 23 cm~ 2 ), obscuring the 
continuum source and allowing the view of only the large distance (hence narrow component) 
of the emission lines. This obscuring structure is referred to as the "AGN molecular torus" , 
considering that it must consist of gas in molecular state, given its low effective temperature 
(T ~ 10 — 100 K). Statistics of Seyfert-1 and Seyfert-2 AGN imply that the height h of these 
torii must be comparable to their distance R from the AGN center, i.e. h/R ~ 1. However, 
the value of this ratio is in conflict with that implied by hydrostatic equilibrium and the ratio 
of their thermal (t> t h ~ 1 (T/100K) km/s) and Keplerian (vk ~ 300 — 500 km/s) velocities, 
namely h/R ~ UthA'K ~ 10~ 3 , thus presenting us with a conundrum concerning the physics 
of these structures. 

These spectroscopically inferred components, along with observations of narrow radio 
jets a long the AGN symmetry axis, led to the now well known AGN picture of ( lUrry fc Padovani 
19951 ). which consists simply of their arrangement at the appropriate positions in the AGN 
vicinity. Compelling as this picture might be observationally, it includes very little, if any, of 
the underlying physics. The AGN constituent components are independent of each other with 
physical properties assigned as needed by the observations of the specific objects. However, 
more recent observational dey e lopme nts suggest that such a picture is rather incomplete. To 
begin with, iBoroson Greenl ( 119921 ) have shown the existence of interrelations among AGN 
the line properties and also relations to other bands of the spectrum (notably the X-rays). 
Then, the increase in UV spectral resolution afforded by HST has shown that roughly 50% 
of Seyfert-ls exhibit UV absorption troughs due to plasma outflowing at v ~ 300 — 100 



km/s , too narrow to have been discerned by the earlier IUE observations ( jCrenshaw et al. 



19991 ). which did detect some, but in a much smaller fraction of the overall AGN population. 
To these flows one must also include those of the so-called BAL QSOs, which reach velocities 
along the observer's LoS in excess of 10 km/s ( iWeymann et al.lll991l ). These are observed 
in about ~ 10% of high luminosity quasars, implying that they subtend a similar fraction of 
the continuum source solid angle in these objects. 

In addition to these UV absorption features, outflowing components were also found 
in the AGN X-ray spectra. The increase in spectral resolution provided by ASCA showed 
that approxi mately 50% of Seyf ert-ls exhibit also blue-shifted absorption features in their X- 
ray spectra ( jGeorge et al.ll2000l ). indicative of outflowing plasma, but of diffe rent ionization 
state th an that responsible for the UV absorption features. More recently, iTombesi et al. 
( I2010af ) have shown that Fe-K absorption features at velocities v ~ 0.1c are rather common 



in nearby Seyfert galaxies and coined for them the term ultra-fast outflows (UFO). The 
simultaneous presence of both UV and X-ray absorbers in the same objects implies they 
belong to the same outflowing plasma (see e.g. iGabel et al.ll2003l ). However, despite a large 
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number of studies supporting this hypothesis, (IMathur et al 



19951 : ICollinge et al.ll200ll : ICrenshaw et al.ll20Q3t Brandt et al 



1994 ; Mathur. Elvis fc Wilkes 



20091 ). an understanding of the 



underlying gas dynamics is lacking. A common origin for the plasma of these components as 
features of a common, radiatively driven flow would be hard to reconcile with their different 
velocities and ionization properties. 

An account of the observed AGN outflows, in particular of the most challeng ing high 
velocity ones of BAL QSOs was put forward semi-analytically by lMurray et al.l ( 119951 ). These 
authors, in analogy with the winds of O-stars, proposed that they are driven off the inner 
regions of the QSO accretion disks by UV and optical line radiation pressure to achieve 
velocities consistent with those observed. The same issue was taken up in more detail in 2D 
numerical calculations by lProga. Stone, fc Kallmanl (120001 ) who included in these calculations 
the detailed photoionization of the line driven wind by the QSO X-ray radiation. As shown 
in this work, efficient wind driving by line pressure requires that the line driven wind material 
be shielded from the ionizing effects of the X-rays, otherwise line driving becomes ineffective. 
Their calculations showed that the "failed wind" from the highly ionized innermost regions 
of the AGN accretion disk did provide the required shielding. The fact that BAL QSOs are 
weak X-ray emitters appears to advocate for such a point of view. 

The ubi quity of AGN outfl o ws im plies that they should be included in the AGN structure 



schematic of lUrry fc Padovanil (119951 ); however, the broad range of observed velocities and 
their different values in the UV and X-ray bands make such a construct complicated in 
the a bsence s of an underlying unifying principle. However, such attempts have been made. 
Thus, lElvisI ( 1200CH ) . motivated by the velocity fields produced by iProga. Stone. &: Kallman 
(120001 ) in model ing the BAL QSO outf lows, proposed a scheme that would supplement the 
AGN picture of lUrry fc Padovanil (119951 ) with outflow components consistent with observed 
phenomenology. By limiting the fast (v <^ 10 4 km/s), radiatively driven flow to a narrow 
angular sliver (A9 ~ 6°) around 9 S ~ 50° (IProga. Stone, fc Kallmanll2000l ). he accounted for 
the observed fraction of BAL QSOs in the overall QSO population. He then attributed the 
lower velocities of the typical X-ray and UV absorption features to the projection effects of 
viewing this flow at a larger angle (8 > 9 S ) and the absence of absorption features in fraction 
of the objects to the low column and high ionization of the wind at 9 < 9 S . He also postulated 
that the angular position 9 S and the opening angle A9 of the high velocity radiatively driven 
radial stream would vary with source luminosity in a way that could account the variation 
of source properties with luminosity. However, his approach ignored the outflows seen in the 
AGN X-ray spectra, which at the time were not as well documented. 



The AGN picture proposed herein is in the same spirit as that of lElvisI ( 120001 ) in that, 
employing a well defined wind dynamical model, it provides a framework of systematizing 
the multitude of observational facts, in particular the more recent high resolution X-ray 
spectroscopy observations. However, it doe s more than tha t; it provides, in addition, scaling 
arguments similar to those put forward by iBorosonl (120021 ) . which allow one to incorporate 
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within this single framework the ionization properties of Seyferts, BAL QSOs and XRBs, 
the structure of the AGN molecular torii and the corresponding IR spectra. This is possible 
because the underlying dynamical models, which span many decades in radius, are to a large 
extent independent of the mass of the accreting object and, as such, they can be applied 
to objects over a very wide range of luminosity. As we will discuss in the ensuing sections, 
the models we present provide the possibility of a broader classification of the structure of 
accreting sources in terms of a small number of parameters (2) thus providing an opportunity 
for a unified treatment of all accretion powered sources. The present work will concentrate 
on the structure of AGN, however it will be argued that the structures of XRBs are quite 
similar, their appearance being different only because of their very different ionization. 

In §2 we provide a brief review of AGN outflow phenomenology with emphasis on the 
more recent high spectral resolution observations of X-ray absorbers by Chandra and XMM- 
Newton. In §3 we present our model and its general scalings. In §4 the model is applied 
to produce the absorber properties of galactic and extragalactic objects as specific cases of 
its parameters along with its general structure properties, depicted in two diagrams that 
relate the absorber column, velocity and observation angle. In §5 we focus on the emission 
properties of these winds and provide an account of the observed linear relation between the 
Ha; and bolometric AGN luminosities as well as their IR-to-far-IR spectra. Finally in §6 
the results are reviewed and directions for future research are outlined. 



2. X-Ray Spectroscopy and Warm Absorbers 



Following the discovery of the ubiquitous nature of X-ray absorption features (Warm 
Absorbers) in the ASCA spectra, the launch of Chandra and XMM-Newton ushered a new 
era in the study of these features. Their superior sensitivity and resolution compared 
wit h those of previous missions made clear the presence of a plethora of transitions (see 
e.g. iBehar et al.l 120031 1 in the spectra of numerous AGN. For example, the X-ray bright 
QSO IRAS 13349+2438 (z = 0.10764) and the Se yfert-1 galaxy (e.g. MCG-6-30- 15 at 
z = 0. 007749) observed e xtensively with ROSAT OBrandt. Fabian fc Pounds! Il996l ) and 



ASCA (IBrandt et al.l 119971 ) gave indications of absorption features in the sub-keV regime. 
Their Chandra and XMM-Newton spectra confirmed this fact and exhibited a wealth of 
transitions such as Nvn, Nevil-Nex, Mgv-Mgxil, Naxi, Sivi-Sixiv, and Nixix while 



successfully identify ing almost all charge states of Fe and O in some cases (IBehar et al.ll2003 



Holczer et al.l 120071 . hereafter HBK07). Furthermore, these observations showed the various 
transitions to be blueshifted relative to the host galaxy, indicating that the warm absorber 
plasma is, in fact, outflowing. 



Though the majority of the X-ray absorber features in the Chandra and XMM-Newton 
archives are associated with the spectra of Seyfert galaxies, absorption features have been 
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also seen in the spectra of quasars; these are mainly Fe-K features usually at velocities higher 
than those seen in Seyferts, in the range of v/c ~ 0.1—0.8, suggesting launching of these winds 
from r egions very close t o the compact o bject . These include b oth BAL [see Chartas et al.l 
d2003h for PG 1115+080. Ichartas et all ( l2007h for H 1413+117. Ichartas et all (120021 l2009h 
for APM 08279+5255; as suc h these also exhibit h igh velocity UV ab sorption featu res 



as 



well as non- BAL quasars [see Pounds et all (120031 ) for PG 0844+349, iReeves et al.l (120031 ) 
for PDS 456, iPounds fc Page! (120061 ) for PG 1211+143]. The columns of these features are 
Nh ~ 10 2 3 — 10 24 cm" 2 which are typically at least ~ 1 times higher than that of UV 
absorbers ( Srianand k, Petit jean 2000 ; Chartas et al.l 20091 ). Finally, as discussed above, Fe- 
K features at velocities v > 10 4 km/ s, have been detected in a large fraction (~ 1/3) of 



observed AGN. iTombesi et al.l (j2010al ) refer to these as ultra- fast outflows (UFO). 



Absorption features o f varying ionization states have also been observed in galactic 
sourc es (GBHC, XRB) (e.g. lMiller et!p200dlNeilsen. Remillard & Ledlioili brandt '& Schulz 
20001 ) . Considering that the presence of the ionic species observed in these sources, like in 
AGN, is due to photoionization of an outflow by the photons of the continuum source, their 
distribution is determined by the photoionization parameter £ = L/nr 2 , where L is the ion- 
izing luminosity, n the local gas density and r the distance of this gas from the source. The 
higher S/N ratios of galactic sources afforded better determination of the velocities of the dif- 
ferent ions as well as the densities of the outflowing wind at a given radius r. Ba sed on these 
param eters, it was concluded (see e.g. §3.1) that the wind of GRO 1655-40 ( IMiller et al. 



20081 ) cannot be driven by the radiation pressure or by the X-ray heating of this source and 



therefore it must be driven by the action of magnetic fields (e.g. iBlandford fc Payndll982 
Contopoulos Lovelace! 119941) . Similar analysis of the X-ray spectra of GRS 1915+105 by 
Neilsen. Remillard fc Leel ~ 2011 ) concluded that the wind mass flux at the outer edge of its 
accretion disk (associated with the lowest ionization ions) could be as much as twenty times 
larger than the mass flux needed to power t he source's X-ray luminosity, while a similar 
conclusion was reached by iBehar et al.l (120031 ) concerning the wind of NGC 3783. Finally, 
it should be noted here that the wind velocities of the highest ionization species (Fe-K) are 
significantly smaller in GBHC than in Seyferts, which in turn are smaller than those of BAL 
QSOs, a fact that should be accounted within the general framework of any wind model. 



2.1. The Absorption Measure Distribution 

The very detection in the AGN X-ray spectra of species of such diverse states of ioniza- 
tion as Fe xxv, Mg v and O I is an important fact in itself. It implies a wind with density 
distribution that produces ionic columns, sufficiently large to be detected, for ions that "live" 
at widely different values of £ and, most likely, also at widely different distances from the ion- 
izing source. This alone is a significant constraint on outflow models. Taking this argument 
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into consideration, iHolczer et al.l ( 120071 ) and subsequently iBeharl (120091 ) fit the ensemble of 
the absorber data in the X-ray spectra of a number of AGN with a continuous distribution 
in £ (rather than adding components of different £ until a sufficiently low \ 2 is achieved). 
In doing so, they developed a statistical measure of the plethora of the transitions in the 
Chandra/XMM-Newton X-ray spectra, the absorption measure distribution (AMD). This is 
the hydrogen equivalent column density of specific ions, N H , per decade of ionization param- 
eter £, as a function of £, i.e. AMD(^) = dN^/d log £. For a monotonic distribution of the 
wind density n(r) with the radius r, determination of AMD is tantamount to determining the 
ionized wind's density dependence on r (along the observers' LoS). Therefore, considering 
that AMD oc Njj, a power law dependence on £ of the form AMD oc Njj oc £ Q , implies also 
a power law dependence for the wind density on r i.e. n(r) oc r~ s with s = (2a + l)/(a + 1). 
So for a _^_0 (AMD independent of £), s ~ 1, consistent with the conclusion reached by 



Beharl (120091 ). 



As noted in IBeharl ((20091), an n(r) oc 1/r density profile is inconsistent with the asymp- 
totic dependence of a mass conserving spherical wind which obtains n(r) oc r~ 2 and whose 
ionization parameter "freezes" to a constant value £oo . On the other hand, considering that 
the velocity of such winds is achie ved only asymptotic ally, with the velocity having the gen- 
eral behavior v(r) = v^l —r^/r) (IMurray et al.lll995l ). the AMD would exhibit in this case 



too a dependence on £ at radii r <^ r*. Assuming a mass conserving spherical wind with the 



above velocity structure, i.e. n(r)v(r)r 2 



m 



dN 



H 



const., one can easily deduce that 

'to 



rn 



1 



dlogi 

where £oo is the ionization parameter at infinity. Under the above velocity field, near r*, 
the wind base, £ (and also v) approaches zero and the gas column diverges so that AMD oc 
N H ^ (m/voor*) (Coo/0- At lar S e distances, with £ = -e, (e/^ < 1), iV^ = {m/voon) e, 
so that the absorption measure distribution has the form AMD oc (e/£<x>) cx r*/r, i.e. the 
wind column decreases like r -1 , as expected. The arguments leading to the above scalings 
are certainly overly simple, however they show that for radiatively driven winds, the column 
decreases with increasing ionization and increa sing velocity , in sig nificant disagreement with 
the dependence found by IHolczer et al.l (120071 ) and IBeharl (120091 ) ; they also suggest that in 
attempting to account for the AMD behavior, one should search for wind models with density 
and velocity profiles asymptotically different from those driven by radiation pressure. 

While it is hard to produce radiation pressure driven winds with the density profile 
implied by the AMD, i.e. n(r) oc r~ l (even an extended source of radiation asymptotically 
appears point-like and the winds acquire density n(r) oc r~ 2 ), MHD winds off accretion disks 



were known to provide such density profiles for some time now (jContopoulos fc Lovelace 
1994 . her eafter, CL94). These are ge neralizations of the two dimensional MHD winds enun- 
ciated by iBlandford Paynd (119821 ) in that they allow also for mass flux rates which de- 
pends on the radius. Given that the velocities of these winds asymptotically scale like 
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the disk Keplerian velocity at the point of launching, i.e. v oc r ' , one can then easily 
see that for n(r) oc r _1 , m(r) oc r 2 n{r)v(r) oc r 1//2 , i.e. the wind mass flux fa increases 
with distance. In this respe ct these winds follow the ADIOS prescription discussed by 
Blandford fc Begelmanl (Il999h . The specific density profile inferred from the AMD is in- 



teresting in that it provides equal column per decade of radius. As such, it is consistent 
with the density needed to account for the lags in the power spectra of XRB and AGN 
(iKazanas et al.l 119971 : iPapadakis et al.l 120011 ) in terms of Compton scattering. Finally, this 
same profile for t he dust density of AGN wind s, at radii beyond t he du st sublimation radius, 
was employed by iKonigl &: Kartjd (119941 ) and iRowan- Robinson! (119951 ) to produce Seyfert-2 
and QSO IR-spectra in agreement with observations. The fact that this wind density profile 
can accommodate in a natural way several diverse and independent aspects of AGN phe- 
nomenology under the same framework is suggestive of its central role in the AGN structure 
and appearance. 



3. The Model MHD Winds 

Motivated by the straightforward interpretation of the AMD in terms of density profiles 
associated with MHD winds, we outline in this section some properties of these wind models 
with emphasis on their scalings with the black hole mass, mass flux rate and distance. We 
pay particular attention on the scalings of the winds' ionization and columns, as these are 
related to the observations of X-ray absorbers. 

The MHD winds considered here are launched by accretion disks threaded by poloidal 
magnetic fie lds under the comb i ned a ction of rotation, gravity and magnetic stresses as 



discussed in iBlandford fc Paynd (119821 ) . To simplify the treatment it is assumed that the 
wind is axisymmetric and as such one need only solve the field geometry and fluid flow in 
the poloidal (r, 6 l )-plane. To simplify further the problem one looks for solutions separable 
in r and 9 with power law r— dependence for all magnetic field and fluid variables. After 
consideration of all conserved quantities in the poloidal plane (mass flux per magnetic field 
flux, angular momentum, magnetic line angular velocity and Bernoulli integral) one is left 
with the force balance in the direction (the so-called Grad-Safranov equation). The 
solution of this equation provides the angular dependence of all the fluid and magnetic field 
variables with their initial values given on the surface of the disk at (r G , 90°), where r Q is a 
fiducial radius that sets the radial scale of the system and it is of order of the Schwarzschild 
radius, r$- The Grad-Safranov equation has the form of a wind equation with several critical 
points, the most important of which for our purposes the Alfven point; this can be crossed by 
appropriately choosing the flow boundary conditions on the disk surface (see CL94, BP82). 

The scalings of the magnetic field, velocity, pressure and density are given by 

B(r,0) = (r/r )-^ 2 B(6)B , (2) 
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v(r, 9) 
p(r, 9) 
n(r, 9) 



= {r/r )-^V{6)Bl , 
= {r/r )- s h{6)Blv- 2 . 



(3) 
(4) 
(5) 



The density normalization is given in terms of the poloidal field on the disk B Q , however, it 
is more instructive to express it in terms of the accretion or wind outflow rate rh as discussed 
in FKCB10; then the density normalization at the inner edge of the disk at 9 = 90° is given 
on setting n(90°) = 1 by 

n o = t> > ( 6 ) 

2a T r s 

where fw is the ratio of of the mass-outflow rate in the wind to the mass-accretion rate 
m , assumed here to be fw — 1, and <tt is the Thomson cross-section. It is important to 
note here that because the mass flux in these winds depends in general on the radius, the 
normalized parameter used throughout this work, m Q , always refers to the mass flux at the 
innermost value of the flow radius, i.e. at r/r = x ~ 1. With the scalings given above we 
then have m(x) = m n x~ s+3 ^ 2 , so for s = 1, as inferred from the fits of HBK07, m oc x 1 ^ 2 



(the iBlandford fc Paynd ( 19821 ) solution has s = 3/2, so n{r) oc x 3 / 2 and therefore m oc x c 
or dm/dlogx = const.). 




20 
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Fig. 1. — (a) Wind density as a function of the polar angle 9 normalized to its value on the disk 
surface (9 = 90°). (b) The line of sight wind velocity on a given field line, normalized to the 
Keplerian velocity of the foot point of the specific field line v Q . 



3.1. Scaling Laws 

The reason that the winds we present here can be applied to AGN as well as to GBHC, 
as mentioned above, is that winds and flows in general are known to be self-similar when: 
The radius r is normalized to the Schwarzschild radius (x = r/rs, r$ = 3M km, where 
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M = M/M is the mass of the accreting object M in units of the solar mass M ); the 
mass flux M is expressed in units of the Eddington accretion rate Me = Le/c 2 (Le is the 
Eddington luminosity, Le — 1.3 ■ 10 38 M erg s _1 ) as rh = M/Me oc M/M; their velocities 
are Keplerian, i.e. v ~ x~ x l 2 c. 

While the wind scalings of eqs. (J2]) - (E]) were introduced for the solution of the MHD 
equations one can write quite general scaling laws that incorporate those of the above equa- 
tions as special cases. We do so now both for accretion flows and winds and use them to 
produce scalings for expressions of the wind density, column density, ionization parameter 
and AMD as a function of the dimensionless radius x, the dimensionless mass flux rate rh 
and black hole mass M. 

1. Accretion: In applying the scalings of eqs. (J2D - © to spherical accretion one 
obtains for the density n(x) oc x~ 3//2 mM _1 and for the column density Nh(x) oc x~ l l 2 rh, 
i.e. an expression independent of M; so, flows onto objects of widely different masses but 
of the same rh = M/Me have the same column at the same distance x from the accreting 
object; its normalization is such to yield Thomson depth tt — 1 for x ~ 1 at m ~ 1 (this is 



the similarity of Advection Dominated Accretion Flows (ADAF) of ( INarayan &: Yilll994j )). 



2. Winds: Analogous considerations of similarity apply equally well to winds, given 
that, generally, their asymptotic velocity, Voo, is proportional to the Keplerian velocity at 
the wind launching radius, r z , i.e. oc x t 1 ^ 2 , with xi = ri/r s . So mass conservation in 
physical and dimensionless units reads respectively 



M ~ n R 2 v 



'oc 



and rhw M oc n(x) x 2 M 2 x l (7) 



with rhw denoting in this case the mass outflow rate in units of the Eddington rate Me- 
Therefore, the wind density and local column density are given in dimensionless units cor- 
respondingly by the expressions (dropping the subscript W from rhw) 

1/2 1/2 

n ( x ) tt-V an d N H (x)(xm— — . (8) 
M x l x 

with the object mass M again dropping out of the expression for the column density Nh- 
For ID winds, e.g. stellar winds or winds driven by the intense radiation pressure near a 
compact object, xi is roughly constant and N H (x) oc rh/x, i.e. the wind column decreases 
inversely proportionally to the distance. However, in the more general case of 2D winds 
(e.g. those of BP82 and CL94) where the wind is launched from a wide range of radii in an 
accretion disk with, generally, different mass flux rates at each radius, x\ ~ x, rh = rh(x) and 
Nh(x) oc m(x)x~ 1 ^ 2 and numerically Nh(x) ~ 10 2i rh(x)x~ 1 ^ 2 cm -2 (it should be understood 
that n(x), Nh(x) have an additiona l angu lar dependence discussed in the previous section). 
As noted in iBlandford fc Begelmanl (Il999h winds with position dependent rh are the means 



by which 2D accretion flows dispose of the excess energy and angular momentum transferred 
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mechanically by the viscous stresses from the smaller to the larger radii, so that finally only 
a small fraction of the available mass accretes on the compact object (BB99). For wind 
density of the form n(x) oc x~ s , the mass flow rate varies with x like (BB99) 

m w (x) = m Q x~ s+3/2 , s < 3/2, (9) 

with m being the mass outflow from the smallest radius x ~ 1.5 — 3; for s = 1, rh oc r 1 / 2 
and N H constant per decade of r. 

3. Photoionization: The wind ionization is determined by the local ratio of photons 
to electrons, the ionization parameter £ = L/n(r)r 2 , (L is the source's ionizing luminosity - 
different from the total luminosity, n(r) the local density and r the distance from the ionizing 
source). This can also be expressed in dimensionless units: If rj (~ 10%) is the radiative 
efficiency of the accretion process, then the luminosity L can be written as L oc i] fn a M 
(m a is the accretion rate that reaches the compact object to produce the luminoisty L), 



or L oc rjfn 2 a M in the case of ADAF (INarayan &: Yil Il994l ) [i.e for m a <, a 2 with a the 



disk viscosity parameter], yielding the for £ an expression also independent of M, implying 
similarity in wind ionization, whether in AGN or XRB (FKCB10) 

Vma - 10 8 -j- — ^— for m a > a 2 (non-ADAF) 



. , - i N H (x)x f w x 3+2 



n(r)r 2 ym a _ 8 y 



N H {x)x fw x 



10 V-zr^ for ™a < a (ADAF) 



where fw = ^o/^ia(~ 1) is the ratio of mass flux in wind and accretion at the smallest 
radii and s the index of the density dependence on r. These relations are also independent 
of M, implying a similarity in the ionization structure of winds, whether in AGN of any 
type or GBHC. For a spherical wind, (s = 2) the ionization parameter is independent of x 
(asymptotically), while for s — 1 it decreases linearly with distance from the ionizing source. 

4- The AMD: Writing Eq. f lTU]) as Nh(x) oc i]m a /l;(x)x one can form the expression 
for AMD (HBK07), namely 

AMD = f H ^\ * 4% (n) 
Gslog£(x) £{x)x 

The fact that AMD is largely independent of £(x) (HBK07) implies £(x) oc 1/x or Njj{x) oc 
log(x), n(x) oc 1/x and m oc x 1//2 , i.e. the wind mass flux increases with radius, as discussed 
in BB99. In this case, both the ionization parameter and the wind density decrease like 1/r 
while the co l umn o f the ions found at a given £ remains roughly constant, in broad agreement 



with iBeharl (120091 ) . So, in disks of sufficiently large radial extent, the wind launched at 
their larger radii will be of sufficiently low ionization and temperature to conform with the 
properties of the AGN unification torus. Finally, one should bear in mind that the above 
scaling laws do not include the 9— dependence of the density and the column density (see 
fig. [T]). At a given (radial) distance r from the source, this leads to much higher ionization 
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of the plasma flowing along the symmetry axis and even smaller neutral absorption column 
than along the higher column density, less ionized eq uatorial directions, in agreement with 



the unification considerations (IKonigl &: Kartjdll994j ). 



Under the above scaling laws, the ionization, velocity and AMD of the winds discussed 
in this work are independent of the mass of the accreting object. Therefore these expressions 
should hold equally well for AGN and GBHC or XRB in general, a fact in obvious disagree- 
ment with observation, given their very different X-ray absorber properties. This would 
indeed be the case if the ionizing spectrum were also independent of the mass of the accret- 
ing object. However, the similarity of wind ionization is finally broken by the dependence of 
the ionizing (X-ray) flux on the mass of the accreting object and leads to the different wind 
ionization properties not only between XRB and AGN, but also amongst AGN: Both in AGN 
and XRB a major fraction of their bolometric luminosity is emitted in quasithermal features 
referred to as the Big Blue Bump (BBB) in AGN and the Multi Color Disk (MCD) in XRB. 
Both are thought to represent emission of the accretion luminosity in black body form by a 
geometrically thin optically thick disk of size a few Schwarzschild radii. The characteristic 
temperature of this feature scales as M _1//4 , being in the X-ray band for solar mass objects 
and in the UV for AGN of M ~ 10 8 M Q . This dominance of the ionizing luminosity (X-rays) 
as a fraction of the bolometric one in XRB over AGN is, according to our model, the main 
cause for their differences in their wind ionization properties as it will be discussed in detail 
in the next section. Even within AGN, the ionizing flux is in fact not a constant fraction of 
the bolometric luminosity but it depends on the absolute source luminosity. A measure of 
their ionizing luminosity is given by the index aox, the logarithmic slope of the flux between 



2500 A and 2 I keV. This quantity has a strong dependence on L(2500 A) ( jSteffen et al.ll2006 



Brandt et al.ll2009l ). a fact that bears relation to the corresponding AGN absorber properties, 



as it will be discussed in the next section. 



4. The Ionization Structure of MHD Winds 

With the wind structure set by the solution of the Grad-Safranov e quation, one can 



now pr oceed to study its ionization structure. The details are given in 



Fukumura et al. 



(j2010al ) so we simply restrict ourselves to a broad description of this procedure. Given the 
(approximately) 1/r density profile, we have split the radius along a given LoS from r D to 
10 6 — 10 8 r Q in approximately 40 segments, spaced equally in logr, so that each has the same 
column density. Then we assume a value for rh which provides the normalization of the 
wind density through Eq. (jfJJ) (we generally assume fw = 1) and a luminosity L related to 
m Q either as L = 1.3 • 10 38 m M erg/s or L = 1.3 • 10 38 m 2 M erg/s (for m < a 2 , according to 
the ADAF prescription) with a specific spectrum; then we invoke the photoionization code 
XSTAR to compute the ionization, opacity and emissivity of the first zone. Using the derived 
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Fig. 2. — The ionization structure of iron appropriate for a Seyfert galaxy with rh ~ 0.1 and an 
ionizing photon spectrum with dF v /du oc z^~ r , V = 1.9, as seen at two different inclination angles 
9 = 30° (left) and 9 = 60° (right). The ionization decreases with increasing distance r along a given 
LoS and with increasing 9 for a given value of r. The column of each ion is given in the left ordinate, 
while the velocity at each position along the LoS in the right one. The LoS velocity at each £ (lower 
abscissa) and the corresponding radius (upper abscissa) is given by the black triangles. 



opacities and emissivity we compute the spectrum exiting this zone, which is used as the 
input for computing the ionization opacity and emissivity of the next zone. The procedure 
is repeated to compute the ionization of the wind along a given line of sight and then for 
different values of the polar angle 9, to produce the ionization structure of the wind over the 
entire poloidal plane. 



Fig. 3 of iFukumura et al.l (l2010al ) depicts the density and ionization structure of such 
winds out to 10 8 r o in linear coordinates. One can see that the wind ionization depends 
both on the distance from the source r and on the angle 9, being highly ionized at small 
values (£ > 10 3 near 9 = 0) even at these large distances, while it decreases to £ < 1 at 
9 ~ 80°. It is also clear that for a given 9, a harder X-ray spectrum will ionize the wind out 
to larger distances. Since the wind velocity along a given LoS depends mainly on the radius 
of launching the specific parcel of gas (but also the inclination angle 9), absorption features 
associated with a given ion should appear at smaller velocities the harder the ionizing spectra. 
This issue is in fact a little more complicated because we only observe the wind velocity along 
our line of sight which depends also on the local shape of the poloidal magnetic field lines. 
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4.1. Ionization of XRB Winds 



The high X-ray content of the spectra of XRB, implies that they should be, in general, 
much more ionized than those of AGN. Nonetheless, X-ray absorption features have also 
been observed in XRBs, as indeed seems to be the case. So P-Cygni Fe-K features were first 
detected in Cir X-l, which is (presumably) an accreting neutron star in a 16.6 day binar y 
orbit around a main sequence companion (IBrandt fc Schula |2000| ; ISchulz fc Brandtl 120021 ) . 
The high ionization of the inner sections of these outflows, render them essentially devoid 
of any absorption features. Therefore, the highest ionization species Fexxvi, Fexxv should 
occur at radii such that £ & 10 4 , which, as the authors suggest sets them at distances r ~ 10 5 
km and corresponding velocities of only v ~ 1, 000 km/s. 

There are several oth er XRB in which simila r features have been detected, the best 

know n being GRO 1655-40 dMiller et al.ll2006l l2008h and GRS 1915+105 flNeilsen. Remillard & Lee 
201 ll ) both of which show Fexxvi and Fexxv absorption features at velocities similar to 
those of Cir X-l. GRO 1655-40 presents a multitude of absorption features similar to those 
of AGN, however, because of the higher S/N ratio their velocities and EW are better deter- 
mined and thus more contraining to the outflow models. Following detailed modeling it was 
concluded flMiller et all l2006l . 120081 ) that the wind of GRO 1655-40 is not consistent with 
being driven by either line radiation pressure or thermally by the X-ray heating of the outer 
regions of the disk. The authors then concluded that the only reasonable means of driving 
the outflow in this object is magnetic, via an MHD wind similar to those outlined in this 
work. 



GRS 1915+105 affords a smilar analysis ( iNeilsen. Remillard &: Ledl201lL and references 
therein). Detailed photoinization modeling of this source provided both the column and 
velocity associated with specific transitions. Assuming that the line absorbing medium has 
thickness Ar ~ r, and that the observed velocities are effectively Keplerian, i.e. v/c ~ 
(V /r,?) -1 / 2 , one can obtain simultaneous estimates of both r and v; these along with the 
measured absorption column Nh — nr yield an estimate of the mass flux of the wind, which 
it was found that it could be as much as several times larger than the mass accretion rate 
needed to power the observed X-ray emission of this source. In fact, one is led to a similar 
co nclusion concerning the mass flux in the wind of Cir X-l by using the parameters provided 



111 



Schulz &: Brandtl (120021 ) . indicating that a wind mass flux increasing with r may not be 



an uncommon feature in accretion powered sources. 

Finally, it is worth noting that the disks of all three of the above sources are apparently 
at high inclinations, i.e. directions that favor observations of such features (less ionization, 
higher column densities). For lower inclination systems, the winds are likely far more highly 
ionized, yielding no apparent absorption features. 
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4.2. Ionization Structure of Seyfert Winds 



The procedure outlined in §4 describing the detailed calcula tion of the ionization o f 
MHD winds was applied to winds associated with Seyfert galaxies ( IFukumura et al.ll2010al ). 
As noted above, given the scale invariance of the wind column and ionization parameter, the 
main difference between the ionization winds off the disks of different classes of accretion 
powered sources (e.g. XRB, Seyferts, BAL QSOs) is the spectrum of the ionizing radiation. 
Since for Seyfert galaxies the X-ray and UV-optical luminosities are approximately the same, 
we have approximated the ionizing spectra with a single power law from 1 — 1, 000 Rydberg 
of photon index T = 1.9 and for the accretion rate we used the value of rh = 0.1. 

The 2D density and p oloidal plane ion i zation structure of a wind for the above param- 
eters is given in Fig. 3 of IFukumura et al.l ( 12010ai ). Here, we present in Figure |2] the ionic 
species distribution of just one element of this wind, namely iron, as a function of the dis- 
tance from the continuum source (upper abscissa) and the corresponding value of £ (lower 
abscissa), for two different values of the observer inclination angle 9 = 30°, 60° (left and 
right panels respectively). The column of each ionic species is given in the left ordinate, 
while the corresponding velocity along the LoS is given in the right ordinate and represented 
by the black triangles. 

We can see that increasing 9 increases the column of a specific ionic species. Also, 
because the radius at which it now forms is smaller, the velocity at the peak abundance of 
this ion is correspondingly higher. Finally, one can obtain the AMD of this flow by adding 
the columns of all ions of Fe at a given value of £. As expected, and shown explicitly in 
(IFukumura et al.ll2010al ). this is constant (i.e. independent of the value of £ or the value of 
r/r$), reflecting the column per decade of radius of the underlying wind. 

The ionization properties of Seyfer t galax ies are consistent with the models presented 
in this figure. The conclusions of iBeharl (120091 ). based on the AMDs of the five AGN he has 
analyzed are broadly consistent with ours as presented in this section. Some of them require 
slightly steeper density profiles (s ~ 1.25) whose ionization structure and velocity properties 
should be similar to those presented in this section. 



4.3. Ionization Structure of BAL QSO Winds 

If the tenets of the proposed model are correct, it should be able to reproduce the velocity 
and ionization structure of the absorption features of accretion powered objects other than 
Seyfert galaxies through judicious choices of the values of the available parameters. Should 
such an enterprise be successful, it could delineate the regions of parameter space occupied by 
individual AGN classes and hence provide their systematization on the basis of an underlying 
physical model. 
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Fig. 3. — (a) 2D poloidal plane density structure of an MHD wind of rh = 0.1 and M = 10 9 Mq, 
ionized by the spectrum shown in (b). The density, log(n[cm~ 3 ]), is shown in color with the 
dotted lines being the iso-density contours, with the log of the corresponding density indicated by 
the associated numbers. The r, z coordinates are in logarithmic scale, so lines of increasing 6 are 
lines at 45° inclination of increasing abscissa intercepts. Also shown are the magnetic field lines 
(black solid) and the loci (white dashed lines) of the positions of formation of Fexxv (circles) and 
Civ (squares), (b) The form of the SED used to produce the wind ionization of (a) consisting of 
(i) a thermal MCD (as BBB) spectrum of innermost temperature kT{ n = 5 eV and (ii) a power-law 
continuum of photon index T. The X-ray flux is normalized relative to that of BBB by a ox . [See 
the electronic edition of the Journal for a color version of this figure.} 



Such an attempt was made by iFukumura et al.l (J2010b|), who employed the specific 
MHD winds to model the absorption features in the spectra of BAL QSOs, in particular 
that of APM 08279+0255 ( jChartas et al"1l2002l . 120031 ). This AGN class was chosen because 
the properties of absorption features of its members are most different from those of Seyferts, 
since they exhibit Fexxv at velocities v ~ 0.5c and Civ at v ~ 0.1c. Because these objects 
were selected on the basis of their broad UV absorption features, any such model should be 
able to provide an account of the combined properties of both their X-ray and UV features. 

The property that most clearly separates Seyferts and BAL QSOs is the relative weak- 
ness of their X-ray to their UV luminosity. This is reflected in the corresponding value of 
their aox parameter, which in BAL QSOs has values apx — —2 compared to — 1 to —1.2 for 
Seyferts. To model the ionization of BAL QSO winds, IFukumura et al.l (j2010bl ) have chosen 
as ionizing spectrum that shown in Fig. |3b- This consists of a BBB component of maxi- 
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mum temperature 5 eV along with a power law component of photon index T, normalized 
at E = 2 keV relative to the BBB so that it results in aox — — 2, consistent with that of 
APM 08279+0255. The decrease in the value of aox is significant, because it implies that 
the ratio of ionizing photon density (given by the X-rays) to the matter density (represented 
by the luminosity of the BBB that peaks in the UV) is quite small, yielding a value for the 
ionization parameter at r ~ r much smaller than that of Seyfert galaxies. Furthermore, 
the cooling effects of the prominent UV emission further contributes to lowering the gas 
temperature at the smallest radii and hence its ionization. It is then possible that some of 
the heavy elements (e.g. Fe) will not be fully ionized even at wind radii r ~ r , leading to a 
high outflow velocity for Fe xxv, Fe xxvi. 

The presence of non-fully ionized Fe at radii r > r a , along with the very steep (un- 
observable) EUV source spectrum, depletes severely the photons at energies E > 50 eV 
needed to produce the ion C +3 of carbon. Since this ion forms at a specific value of £, this 
suppression of ionizing photons is made up by a decrease in the radius r at which it forms 
(and an increase in the local wind velocity), leading to Civ absorption features at velocities 
much higher than those of Seyfert galaxies. Figure [HH shows in logarithmic scale the den- 
sity distribution of a wind of fn = 0.1, and M = 1O 9 M . Dashed lines are the isodensity 
contours with the logarithm of the corresponding density denoted by the assigned number, 
while the thick black lines are the lines of the poloidal magnetic field. The wind flow in the 
poloidal plane is along the magnetic field lines, while the observers' LoS (radial coordinates) 
correspond to lines of 45 degree inclination; an increasing inclination angle 9 corresponds to 
increasing abscissa intercepts for these lines. The white circles and squares along with their 
connecting curves indicate correspondingly the positions of formation of the Fexxv and 
C IV transitions for observers at different values of 9. The ionization structure of the wind 
can thus be translated into the spatial localization of specific ionic species. With a judicious 
choice of 9 (~ 50°) one can obtain r(Fexxv) ~ 5 - 40 r s and r(Clv) ~ 300 - 900 r s (as- 
suming a single LoS) to produce absorption velocities for these transitions consistent those 
of the BAL QSO APM 08279+0255. 

In Fig. H] we show the hydrogen equivalent columns of several ions of iron and carbon 
as a function of the ionization parameter £ for a LoS at an angle 9 = 50° and for the 
ionizing spectrum given in Fig. |3b. This is a figure equivalent to those for Seyferts shown 
in Fig. 2 but appropriate for a BAL QSO. The figure shows also the corresponding wind 
velocity along this LoS (red dashed line), while the shaded regions indicate the positions of 
maximum column for Fe xxv and C IV, which are indicative of the observed velocity of the 
corresponding ion. 
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Fig. 4. — Simulated distribution of local column densities ANh (on the left ordinate) and the 
outflow LoS velocity v/c (on the right ordinate) as a function of ionization parameter £ for (a) C 
and (b) Fe with 8 = 50°. Vertical lines denote the ionization parameter where the local emergent 
column is dominated primarily by Civ in (a) and Fexxv in (b), respectively, also showing the 
corresponding outflow velocity v/c. [See the electronic edition of the Journal for a color version of 
this figure.] 

4.4. The Characteristic Curves of Ionized MHD Winds 

The previous two subsections have shown that, while the underlying velocity structure 
of MHD winds may be very similar and their overall column depending only on the dimen- 
sionless mass flux rh , their ionization structure and the location of specific ions in these 
winds can be quite different depending on their ionizing spectra. Under these (admittedly 
oversimplified) assumptions one can present the kinematic properties of specific ionic species 
of these winds into two diagrams that encapsulate their essential features (namely velocity 
and column) as a function of the ionizing radiation spectral properties and the observers' 
LoS. 

Such diagrams for the transitions Civ(AA1548, 1551A) and Fexxv are shown in Figs. 
[5^.,b. The ClV characteristic curves are shown in blue while those of Fexxv in red. In 
Fig. Et, appropriate for an observer at 9 = 50°, we show: (i) The correlations between 
the ClV, Fexxv velocities and the X-ray photon index F for aox — ~ 2 (solid lines, left 
ordinate), (ii) The correlation between the Civ, Fexxv velocities and the source aox for an 
X-ray photon slope T = 2 (dashed lines, right ordinate). The pairs in parentheses indicate 
the values of (r, aox) at each point. 

We observe that the velocities of both C IV and Fe xxv depend more strongly on aox 
than on T. However, for aox = — 2, as is the case with APM 08279+0255, the Fe xxv velocity 
changes significantly with F and in a way that is actually in agreement with the observations 
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Fig. 5. — (a) The dependence of the model wind LoS velocity (for 9 = 50°) on T (solid lines; left 
ordinate) and aox (dashed lines; right ordinate). The ordered pairs at each point are the values 
of (T,aox) for Civ (blue lines) and Fexxv (red lines). The red squares indicate the change in 
the Fe-K absorption velocity variations with T for APM 08259+5255. (b) The model MHD wind 
correlation between the LoS velocity of a given transition (blue for Civ, red for Fexxv) and the 
corresponding hydrogen equivalent column AiV#(ion) with the observer inclination angle 6 as a 
parameter along these curves (indicated by the numbers in degrees). [See the electronic edition of 
the Journal for a color version of this figure.] 



of IChartas et al.l (120091 ) . given by the squares along the solid red curve. The Civ velocity 
appears to be much less sensitive to T, however it is quite sensitive to the value of aox, 
ranging between v ~ 0.15c and 0.01c as aox varies from —2.1 to —1.6, nicely covering the 
observed range of the velocity of this transition between BAL QSO and Seyferts. 

In Fig. [5b we show the correlation between the LoS velocity of the specific ions (i.e. 
Civ, Fexxv) for a wind ionized by a source with aox = 2, T = 2 (values appropriate 
for a BAL QSO), and the corresponding hydrogen equivalent columns, A?#, with the LoS 
inclination angle as a parameter along these curves. For these specific spectral parameters, 
we see that both velocities increase with 9 up to a critical value different for each transition 
and then level-off as the corresponding ions are produced in parts of the wind that it is 
still accelerating. It is also worth noting that even for the specific spectral distribution, 
appropriate for BAL QSOs, the velocity of C IV is quite small for sufficiently small observer 
inclination angle, while at the same time the corresponding column may be too small for 
the detection of such a feature. Diagrams such as that of Fig. [5b for different ionizing 
spectra and wind parameters will be useful for the determination of the wind parameter 
space (N H ,v) at which specific features become observable. We expect to provide several 
such diagrams for the most important such transitions in future works. 
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5. Wind Emission Properties 

The models and computational setup outlined in the previous were motivated by and 
dealt mainly with the absorption feature properties of accretion powered sources, both galac- 
tic and extragalactic. However, besides the absorber properties, the reprocessing of radiation 
onto these winds provides also emission characteristics, both in continuum and in lines. These 
will be given a first brief treatment in this section, which serves mainly as the road map of 
future work. So in section 5.1 we discuss the effects of dust reprocessing into continuum 
while in 5.2 reprocessing into lines and specifically that of Ha. 



5.1. Torus Structure and Infrared Emission 



The issue of AGN torii and their importance in unification schemes was discussed in 
the introduction. Also outlined there was the conundrum of the disparity between their 
(statistically inferred) height-to-radius ratios (h/R ~ 1) and the values implied by the ratio of 
their thermal to Keplerian velocities (h/R ~ Vth/ yK — 10~ 3 ). The propo sal that resolves this 



conundrum in a straightforward fashion is that of lKonigl fc Kartjd ( 119941 ) who suggested that 
the torii are not structures in hydrostatic equilibrium, but dynamical objects, in particular 
MHD winds driven-off the outer environs of the AGN disks. Furthermore, in fitting the 
observed IR emission as disk UV radiation reprocessed in these winds, they concluded that 
the corresponding density of the reprocessing gas had to be proportional to 1/r as in the 
winds discussed in 83. 



Similarly, iRowan- Robinson! (119951 ) has shown that reprocessing the AGN continuum 
luminosity by dust can provide good fits to their (nearly flat in uF u ) IR spectra, provided 
that the dust density has a profile similar to that used herein, namely proportional to 1/r. 
There is no mystery here: Dust distribution with the specific density profile reprocesses the 
same amount of continuum flux per decade of radius and re-emits it in (roughly) black body 
form; since the dust temperature in regions of increasing radius r falls-off like T oc r _1//2 the 
peak emission will also shift to lower frequencies while emitting the same luminosity; the 
resulting spectrum would th en be "flat" in yF„ u nits, consistent with the Spitzer observations 
of a survey of QSO spectra (INetzer et al.l 120071 ). 



The precise spectrum determination depends, of course, on the details of dust formation 
and properties and also on the details of the reprocessed spectrum and radiative transfer, 
as the emission is not necessarily black body and as the outer layers of the MHD wind 
reprocess not only the continuum luminosity but also the IR emission of the inner "torus" 
regions. One of the more important aspects of the IR and far-IR emission is the frequency 
at which the spectrum turns over from "flat" to the Rayleigh- Jeans form (vF v oc u 3 ), as this 
would provide an estimate of the outer torus radius. Combined space (Herschel) and ground 
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Fig. 6. — Model spectra resulting from reprocessing the BBB AGN continuum (with peak 
at v ~ 10 15 Hz) by winds with three different density profiles s = 0.9 (black), s = 1.0 (red), 
s = 1.1 (blue), that extend roughly two orders of magnitude beyond the dust sublimation 
radius. [See the electronic edition of the Journal for a color version of this figure.] 



(ALMA) observations will be instrumental in this respect. One should simply note here that 
the radial dust temperature dependence, T oc r~ 1//2 , implies that the number of decades in 
radius over which dust reprocessing takes place is twice the number of decades over which 
the AGN vF v spectrum is "flat". In this respect, there should be a synergy between AGN 
X-ray spectroscopy and far-IR observations: The regions at which the far-IR emission turns 
over to the Ray leigh- Jeans form (the "egde" of the torus) are comparable to those at which 
the lowest £ ions should be located; the "edge" of the torus should then manifest then also in 
the X-ray spectra as a steep drop in the X-ray column with decreasing £. Finally, one should 
note that IR spectra similar to those of Seyfert-2 galaxies are not expected in XRB, despite 
their similarity of their winds; this is because the distance at which the dust sublimation 
temperature (T ~ 1000 K) is reached in the wind of an XRB of L — 10 37 erg s _1 , is larger 
than 10 12 cm, the typical XRB binary separation distance. 

Fig. [6] presents an oversimplified example of such a spectrum. In there we present the 
reprocessing of the BBB spectrum of an AGN (the component that peaks near v ~ 10 15 
Hz), by dust which extends from the radius of dust sublimation to the "edge" of the torus, 
a region that spans roughly two decades in radius. The dust density in this region has a 
form n(r) oc r~ s with s = 0.9, 1, 1.1 (black, red and blue curves respectively) to indicate 
the effect that different density profiles on the far-IR AGN spectra. In this specific case the 
AGN is viewed "face-on" and the reprocessing is effectively limited to 9 ^ 50°. The details 
of the X-ray obscuration and its relation to IR emission, bearing relation to the angular 
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distribu tion of our solutions, can be improved with combined X-ray - IR surveys such as 
that of ( jRowan- Robinson et al.l 120091 ). 



5.2. Ha Emission of MHD Winds 



The issue of line emission in AGN is a subject far too broad to be included in any detail 
in the present work. Nonetheless, to our surprise, we found that the basic model we have 
been discussing so far is consistent with well known features and correlations concerning the 
emission line properties of AGN and XRB. So we do make a brief foray into this subject to 
touch upon two issues on which our model seems to provide a novel insight. 

One of the fundamental characteristics of the winds outlined in section 3 are that their 
velocity fields are mainly Keplerian at low latitudes (9 > 70°) while be coming (almost) 
radial at larger ones ( IContopoulos fc Lovelacdll994l ; lFukumura et al.ll2010al ). This particular 
property allows for the possibility of exhibiting emission line features associated both with 
disks and win ds. In fact, there are obs ervations of a "double horn" Ha line profile in 
certain AGN (jEracleous fc Hafpernl 1 19941 ). the tell-tale sign of disk emission. These have 
been modeled as emission excited by the reprocessing of the AGN X-ray continuum on a 
ring-like structure, whose radii and illumination were adjusted to provide agreement with 
observations. The MHD wind models in their most simplistic incarnation, namely that 
described §3, allow for such profiles for the Ha transition, precisely because of the disk-like 
geometry of the flow at low latitudes. Furthermore, the extent of the wind in the 9— direction 
allows also the calculation of the EW of this transition, which is found to be consistent with 
observations (Fukumura & Kazanas in preparation). 



The calculations outlined in §4 provide the photonionization of the entire wind, from 
its smallest (x ^ 1) to its largest (x ~ 10 5 — 10 7 ) scales, set by the size of the accretion disk 
that powers the source. With the wind ionization given, one can also produce scalings and 
correlations involving properties of the corresponding emission of lines of the photoionized 
plasma. While such an in-depth study is beyond the scope of the present paper, some 
straightforward scalings can be produced even at this early stage. With this in mind, we 
compute below the scaling of one of the prominent transitions AGN, namely H a with the 
sourc e luminosity, given that such a relation has already been well documented (lOsterbrok 
19891 ) (we refrain focusing on the Lya and C IV transitions since, being resonant ones, demand 
a more careful treatment within the specific model). 

Quite generally, the luminosity of the Ha transition produced in the MHD winds dis- 
cussed in the previous sections can be written as 

L Ha ~a{T)n 2 Ve Ba (12) 

where a(T) is the recombination coefficient, n 2 V the wind emission measure and eno — 2 
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eV, is the Ha energy. As noted in Eq. (jSJ), the wind density is given by the expression 
n(r,9) ~ n h{Q)x^ l mjM , where h{6) ~ e (e-7r/2)/o.22_ r £\ iB normalization n is set so that 
n r <jT — 1 for m ~ 1 in spherical geometry. For the same value of the total mass flux, since 
the flow in the MHD winds considered herein is concentrated near the equatorial plane, the 
density normalization must be higher by a factor A ~ 1/ f^ 2 h{6) sin# dO ~ 1/0.21 with the 
density profile now reading n(r,6) ~ n a A h(6)x~ 1 rh/M. Then the emission measure of the 
wind reads 

m 2 r n/2 

n 2 V = 2vr n 2 Q — | r 3 M 3 x m 2 / A 2 n(6) 2 sm6 d9 (13) 
M Jo 

~ lOn n 2 r 2 a^ % x m m 2 M ~ 7.1 x 10 49 r o m 2 Mx m (14) 

where the extra factor of 2 in front of the integral takes care of the emission by both hemi- 
spheres, the numerical value of the integral is ~ 2.5, <jt is the Thomson cross section and 
we have taken for the (normalized in units of r$) maximum radius of the disk the value 
x m ~ 10 6 , with r Q ~ rs — 3 x 10 5 cm, the Schwarzschild radius of one solar mass. Then, the 
Ha line luminosity will be 

L Ha ~ a(T) n 2 V e Ha ^ 2 ■ 10 36 m 2 Q M erg s" 1 . (15) 

At the same time, the accretion luminosity will be 

L ~ 10 38 r? m M (16) 

where rj ~ 0.1 is the accretion disk radiative efficiency (77 ~ 0.3 for neutron stars and 
L oc m 2 M in the case of ADAF). 

The obvious point to note between Eqs ( TT5l) and ( TTBl is that they both are proportional 
to the object's mass M; also for m < a 2 (a is the accretion disk viscosity parameter) 
implying that the accretion flow is in the ADAF regime, the ratio R = L^/L of the line to 
the bolometric accretion luminosity is independent of both m and the mass M. This ratio 
depends only on the recombination coefficient a(T) (with assumed value ~ 10~ 13 cm 3 s~ l 
in Eq. ffl5|) above) and the disk efficiency rj and has a value R ~ 1/20 — 1/50. The data 
indicating the correlation between the Ha line and the continuum l umin osity at A4800 A is 
shown in Fig. [71 The data shown are those of 11.6 of lOsterbroki dl989h with the abscissa 



converted from fl ux to luminosity: this figure includes in addition the Ha luminosity of Cir 



X-l (filled square; I Johnston et al.lll999l ). a galactic XRB; inclusion of this point in the figure 
is important because it almost doubles the dynamic range of this relation which now covers 
eight decades in the line luminosity Ln a . The Ha luminosity of this source falls slightly 
below the extrapolation of a linear relation from AGN to XRB. A likely reason for this 
could be the reduced emissivity at small values of x (small radii) du e to the higher gas 



temperature there. Incidentally, for winds with the density profile of the iBlandford fc Payne 
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Fig. 7. — The Ha luminosity of a number of AGN (triangles) as a function of their bolometric 
luminosity (taken to be their monochromatic luminosity at A = 4680 A, as given in O sterbrock 
1989). The square denotes the Ha luminosity of Cir X-l given in I Johnston et al.l (|1999l ) assuming 
the bolometric luminosity to be that of Eddington for a 1.4 solar mass neutron star. 



(119821 ) scaling, i.e. n(x) oc x~ 3 ^ 2 , Eq. ( fl5l) would be proportional to logx m instead of x m 
and the line luminosity would be roughly 10 5 times smaller than observed. We consider the 
straightforward, minimal assumption way that this model accounts for the linear relation 
between the Ha and bolometric luminosities an indication of the validity of its fundamental 
premises. 



6. Summary, Discussion 

In the previous sections we presented a broad strokes picture of the 2D AGN structure, 
one that enc ompasses many decades in radius and frequency and supplements the well known 



schematic of lUrry fc Padovanil (119951 ) with an outflow launched across the entire disk area 
and velocity roughly equal to the local Keplerian velocity at each launch radius. These 
outflows are responsible for the absorption features first detected in the AGN UV spectra 
and more recently established also in their X-ray ones. As noted in the introduction, the goal 
and the spirit of this paper is not to reproduce the detailed phenomenology of specific AGN 
spectral bands, but provide an account of their most general and robust trends; however, 
it purports to do so on the basis of a single, well founded, semi-analytic model for outflows 
off accretion disks that involves a small number of free parameters. In this same spirit, Fig. 
[8] presents a schematic of our model; the radius is in logarithmic space and the shading is 
indicative of the local column density in the spherical— r direction, which is constant, but 
has a strong dependence on the angle 9. The gray lines are illustrative of the magnetic field 
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AGN Disk 



Fig. 8. — A schematic of the model presented in §3. The radius is shown in logarithmic scale, 
with the solid gray lines representing the poloidal field at each radius shown and the corresponding 
wind velocity relative to the fiducial one v Q . The shading is proportional to the local column in 
the r— direction which does not depend on r but has a strong 9— dependence as required by AGN 
unification. At sufficiently large distances, r > 10 4 r a the flow is sufficiently cool to be molecular 
and acts as the torus required by AGN unification. 



geometry with the corresponding scaling of the wind velocity on each one relative to the 
fiducial one v Q shown in the figure. 

As shown in section §3, these outflows have the interesting property that their ioniza- 
tion and dynamical structures scale mainly with one parameter, the dimensionless accretion 
rate m. As such, they are applicable, in principle, to all accretion powered sources from 
galactic accreting black holes to the most luminous quasars. The mass of the object, M, 
simply provides the overall scale of an object's luminosity and siz e (and al s o a ch aracteristic 



temperature for the BBB), in a way similar to that proposed by iBorosonl (120021 ). 



However, despite this economy of parameters, because of the inherently 2D character 
of these winds, their appearance (and that of the AGN central regions) depends quite sig- 
nificantly on the observer's inclination angle 8, a desirable feature and in agreement with 
our notions of AGN unification. Furthermore, and most importantly, as noted in section 
§4, the ionization structure of these flows depends also on the spectrum of the ionizing ra- 
diation. This dependence complicates the situation because it breaks the overall wind flow 
scale invariance on the mass M. As noted earlier, it is effectively the dependence of aox 
on luminosity which is responsible for the differences in the absorption feature properties 
between Seyferts and BAL QSOs. This dependence of aox on L(2500A) (and effectively 
on M) suggests that eventually the wind ionization properties may constitute a two rather 
than three parameter family. Finally, the high ionization of the inner regions of these flows 
in XRB, naturally accounts for the low velocities of the Fe-K features observed in galactic 
sources. 
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The crucial and fundamental aspect of the underlying MHD wind model that allows 
the broad consolidation of the very diverse observational phenomenology of the previous 
sections "under the same roof" is their ability to produce density profiles that decreases like 
~ 1/r with the radius. It is this property that allows the ionization parameter to decrease 
with distance while still providing sufficient column to allow the detection of both high and 
low ionization ions in the AGN X-ray spectra. It is also this property that controls the 
velocities v of the Fe-K transitions in galactic sources, Seyferts and BAL QSO, due to the 
relation between v and £ (v 2 oc £ for the scaling proposed herein; the BP82 scaling leads to 
v oc £ and therefore to much smaller velocity for a given ion, i.e. a given value of £ - and 
also to a much smaller column). Furthermore, it is this specific density profile of the wind 
that allows us to incorporate the physics of the AGN torii within the context of the broader 
physics of accretion onto the compact object, while at the same time producing IR spectra 
in broad agreement with observation. This confluence of the AGN spectral properties with 
the wind spacial structure indicates that AGN and XRB are objects that span many decades 
in radius and frequency, despite the fact that most of their luminosity is released within a 
few Schwarzschild radii. 

The price to pay for a density distribution such as that proposed above is the need to 
invoke winds whose mass flux increases with distance from the source (and more specifically 
like ~ r 1 / 2 ). While this feature can be a ccommodated by the choice of a parameter in the 
models of IContopoulos fc Lovelace! (119941 ). it was given a rather transparent explanation in 
Blandford Begelmanl (1999), in terms of the dynamics of accretion. Interestingly, recent 
Chan dra spectroscopy ( IBrandt Schulzl |2000| ; iBehar et al.l 120031 ; iNeilsen. Remillard fc Lee 
20111 ) appear to support such a notion. At this point, it is not obvious how theoretically 
compelling is this particular feature in the general scheme of accretion properties. Why is 
this mass loss preferred to one that would result in, for instance, the BP82 density profile? 
Are there AGN with winds/accretion flows consistent with n(r) oc r -3 / 2 ? What parameter 
determines which one is chosen by nature? These are pressing questions for which we cur- 
rently have no answers. However, the importance of the specific mass flux dependence on r 
in the interpretation of the AGN X-ray and UV absorber properties will likely attract the 
attention of future studies on this issue. 



Clearly a presentation as broad as that above by necessity ignores a large number 
of issues, both observational and theoretical each of which is in fact a separate branch 
in the study of AGN physics. As such, we have ignored the eff ects of radiation pressure, 
considered in much detail numerically by Pr oga and collaborators (IProga. Stone, fc Kallman 
2000l ; IProgal l2003t IProga fc Kallmanl 120041 ) and semi-analytically by Murray et al.l ( 119951 ) . 
The effects of radiation pressure i n com bination wit h those of the MHD win ds discussed here 
have been considered by lEverettl ( 120051 ) and also by lKonigl fc Kartjd (119941 ) in discussing the 
AGN IR spectra. These works have shown the effects of radiation pressure to be local, as 
they depend on the local flow opacity, thus breaking the similarity of the angular dependence 



-27- 



of the solutions. It was shown in these works that radiation pressure "pushes" to open up 
the field lines to produce a density dependence on 9 different from that given in Fig. [2^ and 
hence would influence the population ratios of objects observed at a given column density. 
However, the effects of radiation pressure do not affect the dependence of wind mass flux 
on the distance (m oc r 1 / 2 ), which is set by conditions on the disk; it is this property that 
determines the radial dependence of the wind density, the property necessary to account for 
their observed phenomenology. 

As noted above, a radiation driven wind, while 2D in the region of launch, it will appear 
radial at sufficiently large distance producing density profiles ~ 1/r 2 there. As appealing 
and compelling as the notion of radiation driven winds is, there is little evidence for them, at 
least in the AGN and XRB X-ray absorber spectra. Advocating radiation driven 2D winds 
across the entire accretion disk, in a fashion similar to those of §3, appears difficult because 
the photon field does not have enough momentum at these large distances to drive a wind 
with the required mass flux. Interestingly the 0— component of the magnetic field, decreasing 
like f?0 cx 1/r has precisely the momentum needed to drive a wind with m oc r 1 / 2 . In this 
respect one should bear in mind that these winds produce most of the kinetic energy at small 
radii, most of the mass flux at large radii and equal momentum per decade of radius. 

Another issue that was only briefly touched upon in section 5, is that of line emission. 
It is generally thought that the line emission in AGN comes from clouds in pressure equi- 
l ibrium with a hot intercloud medium, the result of the X-ray heating thermal instability 



(IKrolik. McKee fc Tarterlll98ll ). Our simple estimates, even though they have ignored this 
possibility, they nonetheless provide an account for the observed correlation of the Ha (a 
transition with minimal radiative transfer nuances) with the AGN bolometric luminosity. 
However, one should bear in mind that our model winds do allow for the formation of such 
clouds at sufficiently small latitudes (below the Alfven surface) where the flow is close to 
hydrostatic equilibrium, i.e. under conditions of a given pressure. We expect that past the 
Alfven point, where the flows are under conditions of given density, the formation of these 
clouds will be less forthcoming. The issue of cloud or wind AGN line emission is an issue 
that deserves more attention and study, given the smoothne ss of the AGN lin es profiles that 



implies a very large number clouds involved in this process ( lArav et al.lll997l ). but certainly 
beyond the scope of the present paper. 

Our treatment of the AGN IR emission has also glossed over much of what it constitutes 
an altogether distinct subfield of AGN study. It has been suggested that clouds are also 
involved in this component of the AGN spectra, however with different properties and at 



radii larger than those of the UV and optical line emitting clouds ( iNenkova et al.ll2008l ). This 
is clearly a point that will have to be looked upon with greater care. A complicating factor 
in this direction is that of star formation in the AGN environment, whose IR contribution 
introduces additional parameters in such a study. Finally, with respect to the IR and far-IR 
AGN spectra, we would like to point to the synergy between X-ray spectroscopy and far-IR 
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observations discussed in §5.1 that would help establish the consistency of this scheme across 
two very different frequency bands. 

Finally, we close with a few words on the "feedback" of our winds on the surrounding 
medium, which likely provides the mass that eventually "feeds" the AGN. First, the angular 
distribution of the flow, as determined by the poloidal field structure, due to its collimation, 
interacts with only part of the surrounding stellar cluster. Second, the energy flux of the 
winds considered here, despite their increasing mass flux with radius, is still dominated by 
the flux at small radii (E = mv 2 /2 oc r^ 1 / 2 ); however, they do carry equal momentum 
per decade of radius (P = mv oc r°), a fact with potential feedback effects, perhaps the 
subject of a future publication. Finally, the radiation effects of the AGN are also limited in 
9 due precise l y to t he winds' column increase with this parameter. However, as shown in 



Tueller et al.l (120081 ) . the fraction of Swift-BAT selected AGN at galactic latitude \b\ > 15° 
with X-ray column Nh > 10 22 cm -2 decreases from ~ 0.5 to close to zero above luminosity 
L ~ 10 44 ergs^ 1 . Therefore, considering that the observed column pertains not only to the 
AGN wind but also to the entire matter distribution along the LoS, the radiative effects 
of high luminosity AGN may affect significantly the evolution of their environment. Such 
constraints will have to be taken into account producing a global evolutionary sequence for 
our models, which are beyond the scope of the present paper. 
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